Anomalous Hall Effect in t 2g Orbital Kagome Lattice due to Non-collinearity: 
Significance of Orbital Aharonov-Bohm Effect 
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A new mechanism of spin structure-driven anomalous Hall effect (AHE) in tilted ferromagnetic 
metals is proposed by taking account of the d-orbital degree of freedom. We find that a conduction 
electron acquires a Berry phase due to the complex d-orbital wavefunction, in the presence of non- 
collinear spin structure and the spin orbit interaction. The AHE driven by this orbital-derived Berry 
phase is much larger than the AHE induced by spin chirality, and it naturally explains the salient 
features of spin structure-driven AHE in pyrochlore Nd2Mo207- Since the proposed AHE can occur 
even for coplanar spin orders (M z — 0), it is expected to emerge in other interesting geometrically 
frustrated systems. 

PACS numbers: 72.10.-d, 72.80.Ga, 72.25.Ba 
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Recently, intrinsic anomalous Hall effect (AHE) at- 
tracts renewed attention as an interesting quantum trans- 
port phenomenon in multiband metals. It is independent 
of the quasiparticle damping rate 7 as shown by Karplus 
and Luttinger (KL) [1]. Recently, the theory of KL had 
been developed intensively, and quantitative studies had 
been performed for /-electron systems [2] , ci-electron sys- 
tems [3, 4], and two-dimensional Rashba electron gas 
model [5, 6]. The origin of giant AHE realized in d(f) 
electron systems is the "orbital Aharonov-Bohm (AB) ef- 
fect" , in which the complex phase factor arises from the 
d(f) angular momentum that is revived by the strong 
atomic spin-orbit interaction (SOI) [2, 4, 7-10]. In d(f) 
electron systems, strong entanglement of orbital and spin 
degree of freedom in multiband Bloch function gives rise 
to the prominent Berry curvature [11, 12]. 

Especially, the AHE due to nontrivial spin structure 
attracts increasing attention, in accordance with the re- 
cent development of the study of frustrated systems. 
The AHE induced by the Berry phase associated with 
spin chirality had been discussed for Mn-oxides [13], py- 
rochlore compounds (kagome lattice) [14, 15], and in spin 
glass systems [16]. Metallic pyrochlore Nd 2 Mo 2 07 would 
be the most famous experimental candidate [17-20]: Be- 
low T c = 93 K, Mo Ad electrons is in the ferromagnetic 
state. Below Tn ~ 30 K, localized Nd 4/ electrons form 
non-coplanar spin-ice magnetic order, and the tilting of 
Mo moment 9 is induced by d-f exchange interaction, 
as shown in Fig. 1. The chirality-driven anomalous Hall 
conductivity (AHC) is proportional to sa • («b x s c) ^ ® 2 
in the weak exchange coupling near half filling [14, 15]. 
However, the neutron diffraction experiments [18] had 
shown that 9 2 <C 10~ 3 , suggesting that the chirality- 
driven AHC is very small. Moreover, the field depen- 
dences of the AHC and 9 2 are quite different [18]. 

In previous studies [13-15], s-orbital models had been 
studied. However, it is natural to expect a significant role 
of the orbital degree of freedom on the spin structure- 
driven AHE, as in the case of the conventional KL-type 



AHE in d-electron systems [2, 4, 7-9]. 

In this paper, we find that the Berry phase is induced 
by the complex d-orbital wavefunction in tilted ferromag- 
netic metals, and it causes a prominent spin structure- 
driven AHE. In Nd 2 Mo 2 7 , the AHE due to orbital 
Berry phase is much larger than the chirality driven AHE, 
since the former is linear in 9 consistently with experi- 
ments. In contrast, both the conventional KL term in 
simple ferromagnets (oc M z oc cos 9) and the spin chiral- 
ity term cannot have ^-linear terms. The present orbital 
mechanism will be realized not only in other pyrochlore 
Pr 2 Ir 2 07 [20], but also in various geometrically frustrated 
metals. 

Here, we construct the i 2g tight-binding model for 
Mo Ad electrons. For this purpose, we consider the 
spinel structure (XMo 2 04) instead of pyrochlore struc- 
ture (X 2 Mo 2 Oy). In both structures, Mo atoms form 
the same pyrochlore lattice. The location of O atoms 
in spinel is rather easy to treat theoretically, and the 
difference is not essential for later study, as we discuss 
later. Figure 1 (a) represents the Mo atoms (black cir- 
cles) and O atoms (white circles) in the spinel structure. 
The [111] Mo layer and the surrounding O atoms in the 
spinel structure is extracted in Fig. 1 (b), where the 
white (gray) circles represent the O atoms above (below) 
the Mo layer. We study this kagome lattice i 2g -orbital 
tight binding model, where the principal axis for the Mo 
d-orbital is fixed by the surrounding Oq octahedron. A 
vector (n x ,n y ,n z ) in the rryz-coordinate in Fig. 1 (a) 
is transformed into [nx ,ny ,nz] m the XKZ-coordinatc 
in Fig. 1 (b) as {n x ,n y ,n z ) — [nx ,ny ,nz]0, where the 
coordinate transform matrix O is given by 
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Arrows in Fig. 1 (c) represents the local effective mag- 
netic (exchange) field at Mo sites, which is composed of 
the ferromagnetic exchange field from Mo 3d electrons 
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FIG. 1: (color online) (a) Spinel crystal structure. Blue 
(white) circles are Mo (O) ions. The Mo ions on the [111] 
plane form the kagome lattice, (b) Kagome lattice made of 
Mo ions. White (gray) circles are O ions above (below) the 
Mo layer. A unit cell contains sites A,B,C. (c) Umbrella like 
local exchange field at Mo sites represented by arrows. 



and the exchange field from Nd 4/ electrons. Under the 
magnetic field parallel to (1,1,1) direction, below the Neel 
temperature of Nd sites, the directions of the local ex- 
change fields at sites A, B and C in the XFZ-coordinatc 
are (<p A = tt/6,0), Ob = 5tt/6,0) and (<fc = 3tt/2,0), 
respectively [17-19]. The tilting angle 8 changes from 
-1.5° (H -► +0 T) to +1.5° (H ~ 6 T) in Nd 2 Mo 2 7 , 
corresponding to the change in the spin-ice state at Nd 
sites [18, 21]. 

Here, we present a discussion of symmetry. If we ro- 
tate the lattice by ir around the site C in Fig. 1 (b), the 
local exchange field is changed from 8 — 8q > to — 8q. 
Since the conductivity tensor is unchanged by the n rota- 
tion, the AHC due to spin chirality mechanism is an even 
function of 8 [14, 15]. However, O sites in Fig. 1 (b) are 
changed by the ir rotation around C, which means that 
the Mo 3d orbital state is changed. For this reason, the 
AHC in pyrochlore compounds due to orbital AB phase 
can have a 0-linear term. 

The i2g-orbital kagome lattice model is given by 

H = ^ t tia,j/3 c ia,<j c j 13,(7 ~~ ^ t h% ' [A t e]cr ) <7'Cj Q , )0 .Cj a!)0 -> 
+A W"./ 3 ' { S U^4a,^,T' ( 2 ) 

where i,j represent the sites, a, j3 represent the t^g- 
orbitals (xy,yz,zx), and a, a 1 — ±1. ti a ,jp is the hop- 



ping integrals between (i,a) and (j,f3). The direct d-d 
hopping integrals are given by the Slater-Koster (SK) 
parameters (ddcr), (ddir) and (ddd). The third term in 
eq. (2) represents the SOI, where A is the coupling con- 
stant, and I and s are the spin and d-orbital operators: 
Their matrix elements are given in Ref. [8]. For conve- 
nience in calculating the AHC, we take the Z-axis for the 
spin quantization axis, where [sx, sy, sz] x 2 becomes the 
Pauli matrix vector. Then, (s x ,s y ,s z ) = [sx,sy,sz}0. 
The second term in eq. (2) represents the Zeeman term, 
where hi is the local exchange field at site i. fi e = —2s 
is the magnetic moment of an electron. 

The Green function is given by a 18 x 18 matrix: 
Gfe(e) = ((e + /i)l — if?.) -1 , where [i is the chemical 
potential and is the matrix for the Hamiltonian in 
the momentum space. According to the linear response 



theory, the AHC is given by ctah = cr AH 
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where G^ A \e) = Gk(e + (— is the retarded (ad- 
vanced) Green function, jk^ = —edHk/dk^ (fi = X,Y) 
is the charge current. Since all the matrix elements in 
jkn are odd with respect to fe, the current vertex correc- 
tion due to local impurities vanishes identically [4, 8] . In 
the band-diagonal representation, eqs. (3) and (4) are 
transformed into eqs. (30) and (32)- (33) in Ref. [4]. 

Before proceeding to the numerical results, we explain 
an intuitive reason why prominent AHE is induced by 
the non-collinearity of the local exchange field hi. For 
this purpose, we assume the strong coupling limit where 
the Zeeman energy is much larger than Fermi energy E-p 
and the SOI. Since fi e = —2s, the SOI term at site i 
is replaced with (— X/2)l ■ n,, where n; = hi/\hi\. Its 
eigenenergies in the t 2g space are and ±A/2, as shown 
in Fig. 2 (a). The eigenstate for E = —A/2 is given by 

1 

-{n x n z - in y )\xy) 



+ (n 2 y + n 2 z )\yz) - (n x n y + in z )\zx)] (5) 

where n = (n x ,n y ,n z ) in the xyz-coordinate: 
(S =A,B,C) in the xyz-coordinate is given by 
[sin 8 cos 4>s, sin0 cos 4>s, cos0]O. When = {n XtVtZ = 
l/v3), eq. (5) becomes \uj\xy) + \yz) + uj*\zx)}/^3, 
where w = exp(i27r/3). Considering that |0| ~ 0.01 in 
Nd2Mo 2 07 [18], we expand eq. (5) up to the first order in 

8. For site C, it is given as | C) = [(I+773 )u;e' l ^ e \xy) + 
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(1 - Jtj.)\yz) + (1 - JL=)u>*e-^ e \zx)]/^. TABLE I 
shows the phases for \a) orbital at site 5 in the eigenstate 
in eq. (5), ip~, up to 0(9). 



TABLE I: Phases for t 2g orbitals in eq. (5) up to 0(9). 
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Here, we explain that a moving electron acquires the 
"orbital AB phase", which gives rise to the prominent 
AHE that is sensitively controlled by the tilting angle 
9. Figure 2 (b) shows the motion of an electron that 
enter into site C via yz orbital and exit via zx orbital. 
When the electron is in the eigenstate |C), the electron 
acquires the phase difference between yz and zx orbitals, 

exp(i(V>°. - ipy Z )) = uj*er l V^ e , which is the orbital AB 
phase that is controlled by 9. The total orbital AB phase 
factor for the triangle path along A^B^C^A in Fig. 2 
(b) is given by the phase of the following amplitude: 

T orb = (A|ff |C)(C|Ho|B)(B|ff |A) (6) 

where Hq is the kinetic term. Here, we take only the 
most largest SK parameter (dda), and neglect (ddir) and 
(dd5) to simplify the discussion. In this case, only the 
following intraorbital hoppings exist: (B; xy\Ho\A; xy) = 
(C; yz\H \B; yz) — (A; zx\H \C; zx) — t. By concentrat- 
ing only on the phase given in TABLE I, we obtain 

T orb ~ ( t 3/ 27 ) e i(^-C)+i(^-0+i(<-^) 

= (t 3 /27)exp(-i3V372(?) +0(9 2 ). (7) 

T or b is also expressed as T or b = |To r b|e -l27r * orb /* , where 
$o = 2irh/\e\ is the flux quantum, and $ or b is the "ficti- 
tious magnetic flux" due to the orbital AB effect. There- 
fore, $ orb = (3/2ir)y^/29 • $0 up to 0(9). Since the 
relation $ rb oc 9 also holds in pyrochlore compounds, 
the orbital AB flux should give rise to the AHE that is 
linear in 9 in Nd2Mo207. The line (i) in Fig. 2 (c) rep- 
resents $ or b for for \9\ < tt/2, which is given by using eq. 
(6) and the relationship -($o/27r)Imln(T orb /|T orb |). At 
9 = ±?t/2, the obtained <I> or b is not an integer multiple 
of $o/2- This fact means that the AHC is finite even in 
the case of coplanar order. 

In the above discussion, we have neglected interorbital 
hoppings integrals for simplicity. If we include them, the 
AHC is finite even if 9 = since the orbital AB phase is 
nonzero, as discussed in Refs. [4, 7]. Thus, interorbital 
hoppings integrals are necessary to realize the AHE for 
9 = (conventional KL type AHE). 

We also discuss the fictitious magnetic flux due to the 
spin chirality, which comes from the Berry phase accom- 
panied with the spin rotation due to the electron hopping 




FIG. 2: (color online) (a) Eigenenergies for t2 g electron under 
the exchange field h and the SOI (-A/2)n ■ I. (b) Orbital AB 
phase given by the complex fe g orbital wavefunction at site 
C. (c) Fictitious magnetic flux due to orbital AB effect (i), 
due to spin chirality (ii), and the total flux (iii). 

[14]. The line (ii) in Fig. 2 (c) represents the fictitious 
magnetic flux due to the spin chirality "I'spin: Geometri- 
cally, — 47r$ sp j n corresponds to the solid angle subtended 
by riA, n-B and no Since ^spin oc 9 2 for 9 <C 1, spin 
chirality mechanism is negligible in Nd2Mo207 [18]. The 
line (iii) shows the total flux <& to t = ^orb + ^spin- 

Thus, all the upward and downward ABC triangles in 
Fig. 1 (b) are penetrated by 4> to t- However, the total 
flux in the unit cell cancels out since the hexagons are 
penetrated by — 2<& to t- Nonetheless, the AHC becomes 
finite [14, 23] since the contribution from the triangle 
path, which is the third order of the hopping integrals, 
would be the most important. The main features of the 
AHE due to orbital AB effect, cr™^, are the following: 
(a) dominates the spin chirality AHE for 9 -c 1. (b) 
er™jj exists even if 9 = tt/2 (coplanar order), whereas spin 
chirality AHE vanishes since $ S pm(# = tt/2) = $o/2. 

From now on, we perform numerical calculation for the 
AHC using eqs. (3) and (4), which contains the effects 
from both <I> or b and 'I'spin, and verify that the concept 
of orbital AB effect is still appropriate even in the weak 
coupling regime where \hi\ <C E-p. Here, we put the 
SK parameters between the nearest neighbor Mo sites as 
(dda) = -1.0, (ddir) = 0.6, and (ddS) = -0.1. | (rfrfcr) | 
corresponds to 2000 K. The spin-orbit coupling constant 
for Mo Ad electron is A = 0.5 [8]. The number of electrons 
per site is n — 2 (1/3-filling) for Nd2Mo207 since the 
valence of Mo ion is 4+. To reproduce the magnetization 
of Mo ion 1.3/^b in Nd 2 Mo 2 7 , we set \h,\ = 0.9. We 
also put the damping rate 7 = 0.2 ~ 0.5 to reproduce 
the resistivity ~ 1 mfkm in a single crystal [19]. 

Figure 3 (a) shows the /i-dependence of the AHC for 
9 = -1.5° (~ T), 9 = (~ 3 T), and 9 = 1.5° (~ 6 
T). The AHC for 9 = is given by the conventional 



4 



n=0 



n=2 



n=6 



i ; » 



= -l.5°(0T) H \ 

!; A 1 

= O°{3T) — //\\ * 



0.5 



<u 

U 
X 

< 



-0,5 




, ('')/ \ 

\ / 

- \ * \ 
\ 
\ 

y '> \ 

£. -i^ 


tb) 

(in) 




(0 

\ 

\ 

\ i 




\ / ' 



0.02 



0.01 



-90 





^[degree] 



90 




1 

& [degree] 



FIG. 3: (color online) (a) AHC for 9 = and ±1.5° as func- 
tion of fi. (b) AHC in the present t2 9 -orbital model for (i) 
7 = 0.5 and (ii) 7 = 0.2. (iii) shows the AHC in the s-orbital 
model for 7 = 0.5. (c) (i) and (iii) for \9\ < 1.5°. 



KL mechanism. Since the present spin structure-driven 
AHE is linear in 9, a very small 9 causes a prominent 
change in the AHC for 1.3 > fi > -1.3 (4.7 > n > 1.0). 
Since e 2 /ha = 10 3 fi -1 cm -1 for a = 4A(a is the interlayer 
spacing), the AHC for n = 2 and 9 — 1.5° corresponds 
to 30 fi cm , which is consistent with experiments. 

Here, we analyze the ^-dependence of the AHC in de- 
tail, by ignoring the experimental constraint \9\ < 1.5°. 
Figure 3 (b) shows the AHCs as functions of 9. The lines 
(i) and (ii) represent the total AHC in the present model 
for 7 = 0.5 and 0.2, respectively. Their overall functional 
forms are approximately odd with respect to 9, and it 
takes a large value even if 6 = ±ir/2 (coplanar order). 
These results are consistent with the AHE induced by 
the orbital AB effect discussed in Fig. 2. Note that 
7 > 0.2 corresponds to "high- resistivity regime" where 
the intrinsic AHC follows an approximate scaling rela- 
tion a ah oc p 2 [2, 4]. The line (iii) in Fig. 3 (b) shows 
the AHC in the s-orbital model studied in Ref. [15]. We 
put the nearest neighbor hopping t = —0.8, 7 = 0.5, 
and \h\ = 0.9. We also put n = 2.5/3 = 0.83 since 
n = 2/3 (1/3- filling) corresponds to band insulating state 
[15]. Although the obtained AHC takes a large value for 
9 ~ 50°, it is 100 times smaller than the AHC in the ti g 
model (i) for |0| < 1.5°, as shown in Fig. 3 (c). 

Here, we make comparison with the theory and ex- 



periments for Nd2Mo2C>7: Under H \\ (1,1,1) below 
In, ct ah decreases with H monotonically from H — > 
T (6 w -1.5°) to 6 T (0 w 1.5°). This experimen- 
tal result is consistent with the AHE due to the orbital 
AB effect that is linear in 9. The present study also 
reproduces the experimental phenomenological equation 

pAHE = 4l] J sM Mo + 47ri ^ M Nd. The firgt term ig 

conventional KL term (9 = 0), and the second term is 



the ^-linear term since 6 oc a/ h 2 x + tiy and yh? x + tiy oc 
M^ d due to noncoplanar magnetic order of Nd moments 
[18, 21]. 

The realization condition for the AHE in the present 
t2 g model is much more general than that in the s-orbital 
model. For example, the AHC is finite even in the copla- 
nar spin state (M z = 0), in high contrast to the chirality- 
driven AHE. Moreover, according to eq. (7), orbital AB 
phase emerges if only one of 9 a, 9b, 9q is nonzero. There- 
fore, in the present ti g model, prominent AHE due to 
non-collincarity is realized unless \\ \\ tiq, with 
the aid of the orbital AB effect and the SOI. 

In summary, we studied the AHE in the ^g-orbital 
model in the presence of non-collinear magnetic configu- 
rations. Thanks to the SOI, local exchange field modifies 
the complex d-orbital wavefunction, and the resultant 
Berry phase induces prominent AHE that is much larger 
than the chirality-driven AHE. Since the derived AHC is 
linear in the tilting angle of Mo moment, 9 , experimental 
large field dependence in Nd2Mo207 is reproduced. 
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